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Abstract 

We study the limiting behaviour of the empirical measure of a system of diffusions in- 
teracting through their ranks when the number of diffusions tends to infinity. We prove that 
the limiting dynamics is given by a McKean-Vlasov evolution equation. Moreover, we show 
that in a wide range of cases the evolution of the cumulative distribution function under the 
limiting dynamics is governed by the generalized porous medium equation with convection. 
The uniqueness theory for the latter is used to establish the uniqueness of solutions of the lim- 
iting McKean-Vlasov equation and the law of large numbers for the corresponding systems 
of interacting diffusions. The implications of the results for rank-based models of capital 
distributions in financial markets are also explained. 



1 Introduction 

The present article studies the behaviour of the weak solutions to the systems of stochastic differ- 
ential equations 

dXi{t) = ii(F j{t) (Xi{t))) dt + <7(F 7(t) (Xi(t))) dBi{t), 1 < i < N (LI) 

on an interval [0, T] in the limit N — > oo. Hereby, j(t) = 4 J2f=i ^Xj(t) is tne empirical mea- 
sure of the particle system Xi(t), . . . , X N (t) at time t, F^m is its cumulative distribution func- 
tion, [i and a are measurable functions on [0, 1] taking values in R and (0, oo), respectively, and 
... , B N are i.i.d. standard Brownian motions. Informally, at any time t the drift and dif- 
fusion coefficients of a fixed particle i are determined by its rank in the particle configuration 
Xi(t), . . . , Xjv(i) at time t, so that whenever a particle changes its rank, the coefficients change 
accordingly. The existence and uniqueness of a weak solution to (11.11) for any N e N was pointed 
out in and is essentially due to the results in @, making this description rigorous. For each 
N E N we fix such a weak solution and denote by Q the probability measure on the space on 
which it is defined. 
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Viewing (ILTT) as an equation for the evolution of the empirical measure j(t) of the particle system 
on [0, T] we prove that under suitable assumptions on /x, a and the initial positions of the particles 
the limiting evolution is governed by a McKean-Vlasov evolution equation. Moreover, in the case 
that the measure in the limiting dynamics is absolutely continuous with respect to the Lebesgue 
measure on E at any time t E [0, T] we show that the cumulative distribution function of the 
system evolves according to the generalized porous medium equation with convection: 

dw d 2 d n . . 

m = 5? sw - T x e{w) (L2) 

where £ and are the antiderivatives of \a 2 and /j, respectively, with E(0) = 0(0) = 0. For 
different values of £ and the equation (lL2l) describes various physical phenomena such as 
infiltration of water into a porous medium or evaporation of water from soil (see E51 and the 
references there). Our law of large numbers for the particle systems in (11.11) shows that the latter 
can be used to obtain numerical approximations of a continuous [0, 1] -valued weak solution of 
the generalized porous medium equation with convection (11.21) provided that E, and the initial 
condition satisfy the assumptions of Theorem 1 .2 below. We conjecture that the same is true un- 
der more general circumstances and, in particular, for any initial condition which is a cumulative 
distribution function of a probability measure. This so-called particle method for numerical ap- 
proximation of solutions of partial differential equations was successfully applied before for the 
Burgers equation (see flU) and the classical McKean-Vlasov equation as in IIT8TI (see B21). 

The system of diffusions in (11711) and related evolutions of particle systems were studied recently 
for fixed values of N 6 N and in some cases for N = oo. They are relevant in the study of capital 
distributions in financial markets (see e.g. 0, lfT9ll . O, [fT3ll and 112210 and in their discrete time 
version in the analysis of the Sherrington-Kirkpatrick model of spin glasses (see e.g. [|20l . 0], 
lETTO . They are also closely related to reflected Brownian motions (see fl~3] for the connection) 
which are widely used as heavy traffic approximations of queueing networks (see e.g. [fTTTl . lfT2l . 

m, USD). 

In the context of capital distributions in financial markets the processes Xi, . . . ,X N stand for 
logarithmic capitalizations of the firms participating in the market. So, our analysis of the limit 
N — > oo of the described systems gives an understanding of the behaviour of the whole market 
under the assumption that the number of firms operating in that market is large. It also allows to 
approximate the evolution of the (logarithmic) capitalization of the j-th ranked firm or of the j 
highest (or lowest) ranked firms in the market for a fixed 1 < j < N under the assumption that 
the number of firms in the market is large. 

Following McKean's seminal work [18J systems of diffusing particles in which the drift and dif- 
fusion coefficients of each particle are functions of the empirical measure of the whole system 
and the position of the particle were studied extensively in the context of particle systems with 
mean field interaction. A good summary of the developments in this direction is given in [|9]|. 
However, in [0 and the references therein the drift and diffusion coefficients are assumed to be 
continuous in the position of the particle and in the empirical measure of the system (with respect 
to the topology of weak convergence on the space of probability measures). This was justified by 
the continuity of potentials and interaction terms appearing in models of statistical mechanics. In 
contrast to this, the coefficients in (11.11) are discontinuous both in the empirical measure j(t) and 
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the position of the particle X^t) which necessitates a more delicate analysis. 



To state the main results of the paper we use the following set of notations. For any separable 
metric space S we denote by Mi(S) the space of probability measures on S endowed with the 
metric 



z: ||z||oo+£«p(z)<l 



d{a, a!) = sup / z(x) a{dx) — I z(x) a'(dx 



S JS 



(1.3) 



where ||^||oo and Lip(z) denote sup^s \z(x)\ and the Lipschitz constant of a function z : S — > 
R, respectively. It is well-known that d metrizes the topology of weak convergence on Mi(S). 
Moreover, we let C([0,T], S) be the space of continuous functions from [0, T] to S, endowed 
with the topology of uniform convergence. We write C([0,T]) for C([0, T], R). For k E N we 
define C*(R) as the space of k times continuously differentiable functions on R having compact 
support, set C£°(R) = P| fc>1 C^(M) and let C C (R) be the space of continuous functions on R 
with compact support, all equipped with the topology of uniform convergence. For any t E 
[0,T] we let Y x {t) < ••• < Y N (t) be the ordered particle system such that Yi(t) = X nt ( i} (t) 
for some (random) permutation n t of {1, . . . , iV} depending on t. In addition, we introduce the 
normalized version of the system X\, . . . , X N given by Ri(t) = X^t) — Y M r N \(0), 1 < i < N 
for t E [0, T] where M(N) = 1 + f if is even and M(N) = ^ if is odd. For any fixed 
iVeNwe let be the distribution of ^ ^jli S Ri(.) as an element of Mi(C([0, T], Mi(R))). 
In the statements on the sequence TZ^ N \ N E N we assume without further mentioning that /i is 
strictly decreasing and choose the initial condition of the particles in such a way that (Y 2 (0) — 
Yi(0), . . . , Vjv(0) — Kjv_i(0)) is distributed according to the unique invariant distribution of the 
process (Y%(t) — Yi(t) , . . . , Y^f(t) — Y/v_i(i)), t E [0,T] (see Proposition 1 of |[T3l for its existence 
and uniqueness). Finally, for each a E M X (R) and / E C C (R) we write (a, f) for J R f da and 
define the operator 

(L a f)(x) = f\x)ti(F a (x)) + l -f"{x)a{F a {x)) 2 (1.4) 
acting on / E C% (R) where F a is the cumulative distribution function of a. 



Our main results can be now stated as follows. 

Theorem 1.1 Let the function fi be continuously differentiable and such that there exists a con- 
stant uq > with fi'(u) < —u for all u E [0, 1] and let the function a 2 be affine. Then the set 
E = {1Z( N ' , N E N} is relatively compact with respect to the topology of weak convergence 
on Mi(C([0, T], Mx(R))) and for any accumulation point TZ°° of "E and any random variable g°° 
distributed according to 71°° it holds 

(Q°°(t),f) - (Q°°(0),f) = I (Q°°(s),L Q ~ (s) f) ds (1.5) 

for all f E C\ (R) and t E [0, T] almost surely. 

Theorem 1.2 Suppose that the functions p, and a are twice continuously differentiable and such 
that for every solution £ of the system 

V/ E C C °°(R) : (m,f) - (Z(0)J) = f (eW.^w/) ds, (1.6) 

m = a (i.7) 
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in C([0, T], Mi(R)) ?/ze measures £(£), t G [0,T] are absolutely continuous with respect to 
the Lebesgue measure on R whenever the initial condition A is absolutely continuous with re- 
spect to the Lebesgue measure on R. TTzen ?/ze solution of the system d/.(5D . d/.7D 15 unique in 
C([0, T], Mi(R))/or eac/j swc/j initial condition. If, in addition, the functions fi and a 2 are affine, 
then the sequence TZ^ N \ N G N converges weakly in Mi(C([0, T], Mi(R))) to to.e Dirac prob- 
ability measure whose atom is the unique solution of the system d/.<5D . d/.7D in C([0, T], Mi(R)) 
wzY/i corresponding initial condition being given in Proposition 3. 1 below. 

Remark 1.3. A careful reading of the proof shows that the statement of Theorem 1.1 is true 
for any sequence of particle systems of the form (1111) provided that [i is strictly decreasing and 
that it holds 




where for any fixed N G N the random vector (Y 2 (0) — Yi(0), . . . , Y N (0) - Y/v-i(0)) is distributed 
according to the unique invariant distribution of the process (Y 2 (t) — Yi(t), . . . , Y N (t) — Kjv-i(t)), 
t G [0, T]. As we show below conditions (lL8l) and (IT9l) are satisfied under the assumptions of 
Theorem 1.1. 

Remark 1.4. In the case that [i is twice continuously differentiable and a is a constant function 
one can deduce the following stochastic representation for an arbitrary solution £ G C( [0, T] , M x (R) ) 
of the system dl.6l) . (H771) whenever the initial condition is absolutely continuous with respect to 
the Lebesgue measure on R. Letting W(t), t G [0, T] be a one-dimensional standard Brownian 
motion on the time interval [0, T] and X(£), t G [0, T] be a weak solution to 

dX(t) = n(F m (X(t))) dt + tr(0) (1.10) 

on [0, T] such that the law of X(0) is given by £(0) one can proceed as in the proof of Theorem 
1.2 below to conclude that for all t G [0, T] the measure £(£) coincides with the law of the random 
variable X(t). Hence, a standard application of Girsanov's Theorem shows that is absolutely 
continuous with respect to the Lebesgue measure on R for all t G [0, T] . Thus, Theorem 1 .2 shows 
that in this case the solution of the system (|L6l) . (11.71) is unique in C([0, T], Mi(R)) whenever the 
initial condition is absolutely continuous with respect to the Lebesgue measure on R. If, in addi- 
tion, the function fj, is affine, then the law of large numbers of Theorem 1 .2 holds for the particle 
systems in (IlTI) . 

Assume now that the functions \i and a are such that the law of large numbers of Theorem 1.2 
applies. In the context of capital distributions in financial markets this means that if the logarith- 
mic capitalizations of the firms in the market follow the dynamics in (11.11) and the number of the 
firms is large, then the evolution of the empirical measure Y^=i $Ri(t) °f normalized logarith- 
mic capitalizations is approximately given by the unique solution of the system (11.61) . (11.71) with 
the initial condition of Proposition 3.1 below. Moreover, the evolution of the capitalization of the 
j-th ranked firm or of the j highest (or lowest) ranked firms in the market can be approximated 
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by the evolution of the jV ~^ +1 quantile or the JV ~^ +1 quantiles (or . . . , quantiles) of 

the solution of the system (ITBT) . (11.71) with the initial condition of Proposition 3.1, respectively. In 
addition, the proof of Theorem 1 .2 shows that the cumulative distribution function of the probabil- 
ity measure 4 ^2f =1 Sr^), which describes the fraction of firms whose capitalizations are below 
a certain threshold, can be approximately described by the unique generalized solution to the 
Cauchy problem for the generalized porous medium equation with convection (11.21) in the sense of 
IfTOll . with the initial condition being given by the cumulative distribution function of the measure 
A°° in Proposition 3.1. 

The rest of the paper is organized as follows. In section 2 we give the proof of Theorem 1.1 
which relies on a characterization of compact subsets of C([0, T], Mi(R)) obtained in J9]|, a char- 
acterization of tight sequences of probability measures on C([0, T]) as in G3l and results on 
convergence of semimartingales in the spirit of [14]. In section 3 we determine the appropriate 
initial condition for the limiting dynamics using Lindeberg's Central Limit Theorem in Proposi- 
tion 3.1 and present the proof of Theorem 1.2 subsequently. In the latter we use a computation 
similar to the one in section 1 of [fT5ll together with the analysis of the Fokker-Planck equation in 
to demonstrate that under the limiting dynamics the cumulative distribution function evolves 
according to the generalized porous medium equation with convection (IT2l) . Using the results of 
IfTOll on the latter we obtain the uniqueness of solutions for the system (11.61) . (|T77l) for any absolutely 
continuous initial condition and as a consequence the law of large numbers of Theorem 1.2. 

2 Proof of Theorem 1.1 

Before proving Theorem 1.1 we recall the results of lfT3l on the invariant distribution of the gap 
process (Y 2 (t) -Y^t),..., Y N (t) - Y N _x(t)), t G [0, T] for a fixed N G N. We remark that the 

results of lfT3l are applicable here, since by subtracting i J2iLi Pyjfjt fr° m Xi(t), . . . , X N (t) 
for all t G [0, T] and by reversing the order of the labels of the particles we can transform the 
particle system in (11.11) into an instance of the particle systems considered in [fT3l . It is shown 
in Proposition 1 of [13] that if /x is strictly decreasing, an invariant distribution exists and is 
unique. Moreover, if the function a 2 is affine, then under the invariant distribution the joint law 
of Y 2 (0) — Yi(0), . . . , Yjv(0) — Fjv-i(O) is that of independent exponential random variables with 
parameters 



Ai(N-i) |2j=i^(/v) N-i X)j=»+i^(iv 



af° = — — '- ^—!— — 1< i < N - 1 

by Proposition 3 of lfT3l . If, in addition, fj, is as in Theorem 1.1, then a straightforward computation 
using the inequalities — A* (if J > — uj for ai l u £ 0, and ^(jA — M n ) > 

shows that 



ujq (u — jjj for all u G 



— 1 



i 5^7=1 ( TV ) jV^Si=i+l' 1 



4 — — > Cn.il) 



a 
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for all 1 < % < N - 1, N G N. 

Proof of Theorem 1.1. 1) By Prokhorov's Theorem it suffices to prove that the sequence 1Z^ N \ 
N G N is tight to show that H is relatively compact. To this end, we fix an arbitrary e > and a 
countable dense subset {/i, / 2 , . . . } of C C (K.) contained in C^(R). From the proof of Lemma 1.3 
in [|9l we see that it is enough to find a compact set K in Mi(R) and compact sets ifi, K 2 , ■ ■ . in 
C([0, T]) such that for all iV G N: 

7^({£ e C([0,T],M 1 (M))|Vt G [0,T] : £(t) G ^ }) > 1 - e, (11.12) 
7^({£ e tf([0,nMi(R))|(£(.),/r) e > 1 - £-2" r , r > 1. (11.13) 

To define we introduce the function (p(x) = a/1 + x 2 and apply Ito's formula to compute 

1 - 

d( e W{t),<p) ^ -^'(R 4 (t))a(F eW{t) (R i (t)))dB i (t) 
i=i 

i=l 

where g^ N \t) = jj X) i=1 ^(t)- The boundedness of </?", fi and cr shows that there exists a 
constant Ci > such that 

1 N 1 

i=l 

for all £ G [0, T] and iVeN. Moreover, for any fixed iV G N the process 
1 - /■* 

Z ^=mJ2 ^(i2 i (s))(T( J P cWW (i2 i ( S )))dB < (s) J £G[0,T] (11.14) 

i=l ^° 

is a continuous martingale. Applying Doob's maximal inequality for non-negative continuous 
submartingales and Jensen's inequality we obtain for all A > 0: 

Q w ( sup (f? (JV) (*),<P) >A + C X T) < QW( sup (6 (Af) + |Z(t)|) > a) 
Me[o,T] ' Me[o,T] ' 

i 

Z(T) 2 5 

where = (£ (JV) (0), <^). Provided that we can show that E Q N [b^] is bounded by a constant 
independent of N, we may employ the Ito isometry to find a constant C 2 > depending only on 
T and sup ue[0)1] er(u) such that for all N G N and A > 0: 

O w ( sup (^(t^j^ + dr) 

v te[o,T] 7 ^ 

Hence, we can choose A such that ^ < £ and let 

K = {a G Mi(R)| (a,<p) < A + dT}. (11.15) 
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As explained in the proof of Lemma 1 .4 in J9J the compactness of the set K in Mi(R) is a con- 
sequence of Prokhorov's Theorem. It remains to show that K Q(N) [b^] is bounded by a constant 
independent of N. It is clear from the definition of £ (JV) (0), N G N that E Q{N) [6 (JV) ] is finite for 
all N G N. Moreover, for TV > 4 we can use (p(x) < 1 + \x\ and inequality (III. 1 II) to compute 

N N 

E« (J >™] = -$>fcp(y<(0) - Fm W (0))] < 1 + ^ $>[|Y,(0) " y W°)ll 
i=i i=i 

Af(JV)-l M(iV)-l JV i-1 

= 1 + AT E E ~W + iV E E ~~w 

i=l j'=i a j i=M(N)+l j=M(N) 3 

9 Sim ( )2 M(N)-1M(N)-1 

- iV ^ 2-r j(iv-j)- 

Finally, the upper bound 

M(N)-1 M(N)-l M(N)-l 

E E 77iv3-y= E F— <MiV-l)-log(iV-M(iV))<log3 

i=l j=i • y '' j=l 7 

for all TV > 4 shows that E Q(JV) [fe^] is bounded by a constant independent of TV. 

2) To prove the existence of sets Ki, PT 2 , . . . with the desired properties it suffices to show that for 
any fixed r G N the sequence of probability measures P^'^, TV G Non C([0, T]) induced by 
1Z( N >, N G N through the mapping £ i— > (£(.), / r ) is tight. To this end, we fix an r G N and aim to 
deduce the tightness of the sequence P^^ r , iV G N from Theorem 1.3.2 of ll23l . To do this we 
need to show 

lim inf P w ' ,r (\y(0)\ < 9) = 1 
etoo NeN 

and 

VA > : limlimsupP (A °' / 'Y sup \y(t) - y(s)\ > A) = 0. 

£4,0 at->oo ^ 0<s<t<T,t-s<e ' 

The first assertion follows immediately by considering 9 > sup xeR |/ r (x)|. To show the second 
assertion we fix a A > 0, define Z(t), t G [0, T] as in step 1, but replacing <p by f r and redefine the 
constant C\ correspondingly. Using the L 2 -version of Doob's maximal inequality for non-negative 
continuous submartingales we obtain for each < e < 

sup \y(t)-y(s)\ > A) < Q {N) ( sup \Z(t) - Z(s)\ > A - de) 



< y: q {n) ( ~p iz( a )-z(i e )i>^-^) 

LfJ-1 _ 2 
^ E (^^J ^ (Ar) [|Z(min((/ + 2),,r))-Z(fe)| 
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where [ J denotes the integer part of a positive real number. The Ito isometry shows that the 
latter expression is bounded by j ( A ~ 2 Cl£ ) 2 ^ with a constant C3 > depending only on 
su Pa;eR \fr(x)\ and sup u£ [ ^ a{u). Taking limits we end up with the second assertion. We con- 
clude that the sequence 1Z^- N \ N E N is tight. 

3) To prove (11.51) we let IZ 00 be the limit of a converging subsequence TZ^ Nk \ k E N of the 
sequence TZ^ N \ N E N. Next, for each k E N we let ff Nk > be a random variable with distribu- 
tion TZS Ni >) and g°° be a random variable with distribution 72.°°, all defined on the same probability 
space and such that g^ Nk > — )^oo £?°° in C([0, T], Mi(R)) almost surely. This is possible due to the 
Skorohod Representation Theorem in the form of Theorem 3.5.1 in Q. Indeed, the metric space 
C([0, T), Mi(R)) is separable, since the countable set of functions whose values at 0, j, 4, . . . , 1 
belong to a fixed countable dense subset of Mi (R) and which interpolate linearly on the intervals 
foraj € Nis dense in C([0, T], Mi(R)). 



0,4 



3 ' 



From the dynamics computed in step 1 we observe that for each k E N and / G C^(R) the 
process (g^ Nk \t), f),t E [0, T] is a semimartingale in the sense of definition II. 2. 6 in llT4ll with its 
characteristics being given by 



f'(x)fi(F^ k)(s) (x)) + -f"(x)a(F^ Nk){s) (x)) 2 # N »( 8 )(dx) ds, 

jr f j f^f<F^ k){s) {x)f # N «\s){dx) ds 

which we denote by B^ Nk \t) and C^ Nk \t), respectively. We claim that in order to establish (|T3 
it suffices to show that for any f E (R) we have 



sup 

te[o,T] 



'3F(s),Lp°( a )f) ds 



Ao, 



E[tfW>(T)] 



(11.16) 
(11.17) 



where — > denotes the convergence in probability. Indeed, the second convergence together with 
the L 2 -version of Doob's maximal inequality for non-negative continuous submartingales would 
imply 



sup 

te[o,T] 



sup 

te[o,T] 



:d iNk \t)j)-^ Nk \o)j)-B^ k \t) 4o. 

Hence, from the first convergence we would be able to conclude that 

h ds 



(11.18) 



(^(0),/) 



4o. 



(11.19) 



By a diagonalization argument relying on the separability of C^(R), endowed with the topology 
of uniform convergence of functions and their first and second derivatives, we would be able to 
find a subsequence of Nk, k E N such that the latter convergence holds for all / E C\ (R) in the 
almost sure sense. But since (]f Nk \t), f) — (g^ Nh \0), /) — f*(g°°(s), L^^f) ds converges to 
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(Q°°(t)J) ~ (q°°(Q)J) - f5((?*'(.s),Lp°(.,)f) ds in the limit k -»■ oo for all / G C C 3 (M) and 
i £ [0, T] almost surely, we would obtain equation (II3T) . 



4) We now show the two claimed convergence results. The convergence in (III. 171) is a direct 
consequence of the boundedness of /' and a. To prove the convergence in (III. 161) we introduce 
for any e > and ieRa Lipschitz function j e x such that 



1(— oo,x— e] — fx — 1(— oo,x] 



(11.20) 



and /J, is a translate of /J by x' — x for any x,x' G E and e > 0. Moreover, for any probability 

measure q on 1 we set F£(x) = {a, f e x ). We note B^\t) = / *(^ fc) ( s )» L ^) {s) f) ds for all 
t G [0, T] and fcGN and deduce from the triangle inequality that 



sup 

te[o,T] 



< sup 

te[o,T] 



+ sup 

te[o,T] 



[ @ Nk Ks),L^ Nk){s) f) ds- f\r(s),L^ {s) f) ds 
Jo Jo 

@ Ni) {s), f{n o Fg( Wjk , w )) - (fiP°(«), f{li o Fy-w)) ds 



where o denotes the composition of functions. We claim that the latter two terms converge to zero 
in probability. Since the proof of this claim is identical for both terms, we only carry it out for the 
first (drift) term. To this end, we fix an e > and a k G N and observe 



sup 

te[o,T 



[ (^(s), /'Gu o F Ws) )) - (F°(s), /'(/x o Fy» w )) ds 
Jo 

f \$ Nk \s\ f(j* o i> fc , w )) - /'(/. o FJ Wfe)(s) ))| ds 

jf (^H*). f'Qi o F^ fc)(s) )) - (g»( a ), f Oi o Ff^)) ds 

</ 



< 



+ sup 

te[o,T] 



We call the summands on the right-hand side (I), (II) and (III). We will bound the three terms 
consecutively. 



Denoting by Y±(t) <••• < Yff (t) the ordered particles of the normalized system R\(t), . . . , i?jv & (t) 
for any time t G [0, T] we can bound term (I) from above by 



r T i Nk i i 



where the constant C4 is the product of sup.j, gR \f'(x) | and sup ug [ ^ I • By Fubini's Theorem 
and the defining property of the initial condition (see the paragraph preceeding Theorem 1.1) it 
follows that the expectation of (I) is bounded above by X£i Ej=i F (^(0) - 1^(0) < e). To 



bound this expression further we first choose a C5 > such that a\ < C$N for all N £ N and 
1< i < N - 1. This is possible due to the obvious bound aJ 7V) < • mp « e ™ which 

— — f i — N inf„ 6[0 ,i] v{uy 

holds for all AT G N and I < i < N — 1. Next, for any fixed k G N we let E u . . . , E Nk _ x be 
i.i.d. exponential random variables with parameter C 5 N k and P be a Poisson random variable with 
parameter C 5 N k e. Then from the scaling property of exponential random variables we deduce 



Ni. 



% £ £ (°) - m < *) < ^ E E w + • ■ ■ + e ) 



m 



C4T 



»=i i=i 

EE'(^h)<|(i + w)^(. + »). 

i=i j=i 



i=i i=i 



C 4 T 

4 



All in all, we conclude that 

r r T 

lim sup E / 

k— >oo \_Jo 

is bounded above by C 4 C 5 Te. 



\(# N "\s), /'(/x o F^ ft)(s) )) - (^( S ), ftp o FJ Wfc)(s) ))| 



Term (II) can be estimated from above by 



sup 

te[o,T] 



+ sup 

te[o,T] 



./o 



Now, recalling the definition of C 4 we bound the first summand from above by 



<e[o,T]Jo 



C 4 sup / sup|F^ )w (x) - F~oo {s) {x)\ ds. 



(11.21) 



Since for each x G K the function /J is Lipschitz with the Lipschitz constant being independent 
of x and the convergence d(jf Nk \s), £°°(s)) — >k->oo is uniform in s, the expression in (III. 2 II) 
tends to zero almost surely in the limit k — > 00. Moreover, for any x, x' G K it holds 



ff- w GOI < sup \f x (y) - f x (y -x + x')\ < Lip(/«) ■ \x - x'\. 



Hence, the uniformity in s of the convergence d(g^ Nk ^(s), — hfe->oo implies that the second 
summand in the bound on term (II) tends to zero almost surely in the limit k — >■ 00. Thus, from 
the Dominated Convergence Theorem we deduce that the expectation of term (II) converges to 
zero in the limit k — > 00. 

Due to the inequality 



- F^ {s) (x)\ < F§ 



e°°(s)' 



x + e 



- F~ 



(11.22) 
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for all x e K the expectation of term (III) can be bounded above by 



C A E 



f I (F§ aa{s) (x + e)-F§ 00{s) (x))r , (s)(dx) ds 
Jo Jr 



Moreover, applying the triangle inequality as in the proof of the upper bound on term (II) one 
shows that with probability 1 the integrand in the cfe-integral is the limit of the corresponding 
objects with replaced by $ Nk \ Indeed, the difference between the two is bounded above in 
absolute value by A(Lip(fo) + 1) • d(g^ Nk \s), g°°(s)). We conclude from Fatou's Lemma that the 
expectation of term (III) is bounded above by 



C 4 lrminfE 

fc— >oo 



T 



(x + e)-Fl iNk){s) (x))Q( Nk \s)(dx) ds 



Bounding the integrand in the rfs-integral in the same way as in the corresponding estimate on 
term (I) one obtains the upper bound 



C 4 lim inf E 

k—too 



ds 



Finally, proceeding as in the upper bound on term (I) we deduce that the expectation of term (III) 
is bounded above by 2CiC 5 Te. 



All in all, we have shown that 



lim sup E 

fc— >oo 



sup 

te[o,T] 



f'(f, o i> fc)(s) )) - (r>(s), /V o F^ {s) ))ds 

Jo 



is bounded above by Cqe with Cq = 3C4C5T. By taking the limit e I we deduce that the 
term inside the latter expectation tends to zero in L 1 in the limit k — > 00 and so, in particular, it 
converges to zero in probability. □ 



3 Proof of Theorem 1.2 

To be ready to prove Theorem 1.2 we show next that the initial probability measure £(0) is the 
same under each accumulation point 1Z°° of the set S provided that the functions \i and a 2 are 
affine. 

Proposition 3.1 Let the functions fi and a 2 be affine, so that <j 2 (u) = cu + dfor some constants 
c,d£R Then for any choice ofTZ°° as in Theorem 1.1 the distribution of the initial probability 
measure £(0) under TZ°° is a Dirac probability measure. Moreover, the atom of the latter is given 
by the unique A°° G Mi(M) whose quantiles q°°(u), u e (0, 1) are given by 

^ ( " ) =-^ log(2 - 2 " )+ ^k log(2i,) - (IIL23) 

In particular, A°° is absolutely continuous with respect to the Lebesgue measure on ffiL 
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Proof. 1) For each N G N let q( N ~> (u),u G (0, 1) be the (random) quantiles of the initial probability 
measure £(0) under TZ^. Then for any fixed \ < u < 1 it holds 



|"uAT|-l 

q {N \u)= (^ + i(0)-^(0)) 

i=M(N) 

_Et% 1 ) (Y, +1 (0)-Ym-rn^ 



for all iV G N such that \uN] - 1 > M ( JV) where 



m 



(JV) 



|"uAT|-l 

E 

i=M(N) 



a 



1 

(AO' 



,(A0 



E 

i=M(N) 



(«T) a 



(111.24) 



and denotes the smallest integer greater or equal to x for any x G R. Using the definition of 
af^, l<i<JV-l,JVeN and the assumption that jj, is affine we see that for any | < -u < 1 

there exist constants wi(-u), w 2 (-u) > such that ui(u)N < a\ N ^ < u 2 {u)N for all M(N) < i < 
\uN] — 1, N G N. From this and Lindeberg's Central Limit Theorem we deduce that for any 
fixed \ < u < 1 it holds 



lim q^ N \u) = lim 



Af-5>oo 



Af-5>oo 



(111.25) 



in distribution, provided that the latter limit exists. Indeed, we can check Lindeberg's condition 

y i+1 (o)-yi(o)-^ 

for the triangular array ^ (N) '—, M(N) < i < \uN~\ — 1, N G N as follows: for any 

e > and any N G N such that \uN] - 1 > M(N) we have the estimates 



iwi-i 



£ E [(^+l( ) - Y i( Q ) ~ -Jn)) ' \Y^0)-Y i (p)--X n \>eVm] 



i=M(N) 



\uN]-l 

i=M(N) y u i ) 



(^ r )^2 



• E 



G 



< E 



G 



M(N),N 1 -\G MWtN \>ew 3 (u)VN 



where we have set G itN = (Y i+1 (0) - E;(0))af° - 1 for M(N) < i < \uN] - 1, N G N 
and w 3 (u) > is a constant such that Vv^af 1 > ou 3 (u)y/N for all M(N) < i < \uN] - 1 
and N G N. The latter expectation tends to zero in the limit N — > oo, since the random vari- 
able Gm(n),n + 1 is distributed according to the exponential distribution with parameter 1 for all 
N > 3. Hence, Lindeberg's condition is satisfied. 
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2) Now, plugging in the definition of aj for M(N) < i < \uN] — 1 and N E N we calculate 

[uN]-l q+c/2 



lim nv- N ' = lim 



N c -^ + d 



iV-Kx> TV-s-oo ■ / — ' i(JV — i) lit' 

i=M(N) K ' ^ 

\vN\-l , \uN]-l 

c 



lim > h ^ — — — r • lim > ( — I ) 

\U'{[})\ N-^oo ^ N-i Lu'(O) iV->oo ^ W N-i 
lr V 71 i=M(N) V 71 i=M(N) V 7 

.. /c + tZ. /N-M(N)\ d , / uN \\ 

J^o 1 17(5)1 log I iv-niv j + Ro)l g ImPO ^ ■ 



Recalling that y < M(N) < | + 1 we can compute the last limit to 

c + d , . „ . d , . . 

~Ro)l g ~ + s 

An analogous application of Lindeberg's Central Limit Theorem and similar calculations to the 
ones above show for the case < u < |: 

M(JV)-1 1 c + rf ^ 

lim o^fn) = — lim > —7^ = —-. — — T log (2 — 2u) + 1 — r-rr log(2u) 

K) ^4^°* H0)l l//(0)l 

in distribution. 



3) If IZ 00 is as in Theorem 1.1, then the Skorohod Representation Theorem in the form of The- 
orem 3.5.1 in Q shows that we can find an increasing sequence N k , k E N of natural numbers 
and random variables ^^(0), k E N and £?°°(0) defined on the same probability space such that 
for each k E N the distribution of the random variable ^^(O) is given by the law of £(0) under 
Tl( Nk) , 0°°(O) is distributed according to the law of £(0) under TZ 00 and g {Nk) (0) -► fc - +00 £°°(0) 
weakly with probability 1. It follows that the quantile functions of ^^(O) converge in the limit 
k — > 00 to the quantile function of g°°(0) at all continuity points of the latter almost surely (see 
e.g. the proof of Theorem 2.2.2 in [8]). By Fubini's Theorem we obtain that the w-quantile of 
gi N k)^ converges to the w-quantile of £?°°(0) almost surely in the limit k — > 00 for Lebesgue 
almost every u E (0, 1) and, in particular, for all u in a countable dense subset of (0, 1). Due to 
the monotonicity of quantile functions and the computations in steps 1 and 2 the quantile function 
of £?°°(0) has to coincide with q°° (defined in the statement of the proposition) with probability 
1. Hence, the distribution of g°°(0), which is the same as the law of £(0) under 1Z°°, is given by 
the Dirac probability measure described in the proposition. Finally, the probability measure A°° 
is absolutely continuous with respect to the Lebesgue measure on M, because its quantile function 
q°° is continuously differentiable and strictly increasing on (0, 1). □ 

Combining the ideas of [|9]| and lfT5ll with a result in IfTOll we can now prove Theorem 1.2. In 
the proof we use the following notations. For a measurable subset S of a Euclidean space we 
write L' P (S) and ||.||lp(s) for the space of functions / : S — > M such that |/| p is integrable with 
respect to the restriction of the Lebesgue measure to S and the corresponding L p -norm, respec- 
tively, where p is a real number in [1, 00). In addition, for any real- valued random variable Y we 
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denote by C(Y) the law of Y. 



Proof of Theorem 1.2. 1) It suffices to prove that the initial value problem (ITBT) . (11.71) has at 
most one solution in C([0, T], Mi(R)) for each initial condition which is absolutely continuous 
with respect to the Lebesgue measure on R, since then the law of large numbers follows from 
Theorem 1.1 and Proposition 3.1. To achieve the former we fix an absolutely continuous initial 
condition A, let v u i G {1, 2} be two solutions of C6]), dO) in C([0, T], Mi(R)) with this initial 
condition and will prove v\ = u 2 . We show first that for i G {1, 2} the measures z^(t), t G [0, T] 
are given by one-dimensional distributions of solutions of appropriate martingale problems. To 
this end, for i G {1, 2} we define Zi(t), t G [0, T] as the respective unique (in law) solutions of 
the martingale problems associated with the families of operators L u .^, t G [0, T] such that for 
any / G C£°(R) the processes 



c 



/(^(t)) - /(Zi(0)) - / (^ iW /)(^i(s)) ete, t G [0,T], 

Jo 

i G {1,2} are martingales and C(Z X (G)) = C(Z 2 (0)) = A. Due to Exercise 7.3.3 in HI the 
processes t G [0, T] and Z 2 {t), t G [0, T] are well-defined and C{Zi{.)) G C([0, T], Mi(R)) 

for i G {1,2}. We claim that z/j(.) = for i G {1,2}. To prove the claim we fix an 

i G {1, 2}. By their respective definitions Ui(t), t G [0, T] and C(Zi(t)), t G [0, T] solve the initial 
value problem for the Fokker-Planck equation 

V/ G C C °°(R) : (£(*),/) - (£(0),/) = f\t(s),L Ma) f) ds, (111.26) 

Jo 

^(0) = A (111.27) 

on [0, T]. We show now that the latter has a unique solution in C([0, T], Mi(R)). To this end, we 
observe that the operator 1Z = ^ + A^q is continuous as an operator from the Sobolev space 
VK 1,2 ' p ([0,T] x [— r, r]) (the space functions in L p ([0,T] x [— r, r]) whose generalized first time 
derivative and generalized first two spatial derivatives belong to L p ([0,T] x [— r, r]), endowed 
with the usual Sobolev norm) into L p ([0, T] x [— r, r]) for any p > 1 and any r > 0. This is due 
to the boundedness of /x and er. Moreover, following the steps in the proof of Theorem A.l in 
we obtain for any p > 3 and r > 0: 

(^(0, + L Mt) )f) dt > (c(Z t {t)), + L Mt) )f) dt (ni.28) 

for all / G W 1,2,p ([0, T] x [-r, r]) (the space of functions in W /1 ' 2,p ([0, T] x [— r, r]) vanishing on 
([0, T] x {— r, r}) U ({T} x [— r, r])) such that 7?./ > Lebesgue almost everywhere. Hereby, we 
have used the convention {^+L Vi {t))f = on the complement of [0, T] x [— r, r] in [0, T] x R. By 
Theorem 9.1 in chapter IV of IfTTl the image of the just described functions under 1Z is given by 
lq_([0, T] x [— r, r]), the set of functions in L p ([0, T] x [— r, r]) which are non-negative Lebesgue 
almost everywhere. Indeed, from Proposition 3.1 and the assumption of the theorem we conclude 
that for every t G [0, T] the measure Ui(t) is absolutely continuous with respect to the Lebesgue 
measure on R. Thus, the function (t, x) (-)■ a(F Vi ( t )(x)) is continuous and Theorem 9.1 in chapter 
IV of IfTTll is applicable. Since r > was arbitrary, we deduce from the Monotone Convergence 
Theorem that 

cT p 

g(x) ■ l[ ,t](s) Ui(s)(dx) ds> I I g(x) ■ l[ 0)t ](s) C(Zi(s))(dx) ds 
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for all non-negative continuous bounded functions g on E and all t G [0, T]. Using the same 
inequality with (sup xen g(x)) — g instead of g we infer that equality must hold in the latter in- 
equality. Differentiating the resulting identity with respect to t we see that Ui(t) = C(Zi(t)) for 
all t G [0,T]. 



2) To finish the proof we aim to apply Theorem 4 of IfTOl . Approximating the functions f\ i 



x 



x 



and f 2 (x) 



x 



by functions in C£°(R) coinciding with fi, f 2 on [—A, A] for increasing values of 



AeN and applying Proposition 4.6 in chapter 5 of [|T6ll we conclude that there exist probability 
spaces on which processes of the same law as Z\ and Z 2 (which we will also denote by Z\ and 
Z 2 ) are defined such that 



dZi(t) = //(F„. (t) (£<(*))) dt + a(F Mt) (Z t (t))) dW { {t), t G [0,T] 



(111.29) 



holds for i G {1,2} and appropriate standard Brownian motions W\, W 2 . Next, we fix arbitrary 
numbers x x < x 2 in E and £ x < £ 2 m [0, T] and introduce the function f(t, x) = ip(t, y) dy 
on [ti, t 2 ) x E where ^ is an arbitrary continuous function on [ti, t 2 \ x E which is continuously 
differentiable in both variables with ?/>(., sc) = = whenever a; ^ (xi,X2). Applying Ito's 

formula to f(t, Zi(t)) and taking the expectation we obtain 



(z/ J (t 2 ),/(t 2 ,.))-(z/ l (t 1 ),/(t 1 ,.)) 
f t2 ( df df 



»)(•)) + 5^(t,-M^(.)0) a )'« 



for z G {1,2}. Recalling that i/i and v 2 are solutions of the Fokker- Planck equations in step 1 and 
following the proof of Lemma A. 2 in [9 ] we conclude that the finite measures corresponding to the 
functionals h H- J Q T L h(t, x) i>i(t)(dx) dt, h i— >■ J Q T J K x) u 2 {t){dx) dt acting on continuous 
bounded functions on [0, T] x E are absolutely continuous with respect to the Lebesgue measure 
on [0,T] x E. Moreover, the proof of Lemma A.2 in |[9l shows that the corresponding density 
functions k%, k 2 on [0,T] x E are locally square integrable. Setting Wi(t,x) = F n (^{x) for 
t G [0, T], x G E, i G {1,2} and applying integration by parts with respect to the spatial variable 
in the last equation (recalling from step 1 that the measures Ui(t), t G [0, T], i G {1,2} are 
absolutely continuous with respect to the Lebesgue measure on E) we see that for i G {1, 2} it 
holds 



ip(t 2 , x)u)i(t 2 , x) dx — 
t2 r /dip 



ip(ti, x)wi(ti,x) dx 
dip 



t2 



^— (t, x)wi(t, x) + — (t, x)Q(wi(t, x))j dx dt 
dip 



dx 



(t, x)H'(wi(t, x))ki(t, x) dx dt 



which we will call equation (*). Hereby, and S are the antiderivatives of // and |cr 2 , respectively, 
for which 6(0) = S(0) = 0. Next, we note that for any n G N and any continuous function 
ip : [ti, t 2 \ xE — >■ E which is supported on D = [ti, t 2 ] x [xi, x 2 \ and is continuously differentiable 
in t and twice continuously differentiable in x on D we can find a function ip n of the same type as 
the function ip in equation (*) such that ip n — >„^oo ip uniformly on D and 



< -. 

l 2 (d) n 







dlp n 


dip 




dlpn 


dip 


+ 


dt 


~ ~dt 


+ 

L 2 (D) 


dx 


dx 
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This can be achieved by modifying ip on small neighborhoods of 

QpiM x {xi}) U {[t u t 2 ] x {x 2 }) 

in D. Thus, an approximation argument together with the Cauchy-Schwarz inequality shows that 
equation (*) holds for all functions ij) of the described type and any i E {1, 2}. Choosing fc, such 
that for each t E [0, T] the function ki(t, .) is a density function of the probability measure u^t) 
for z e {1, 2} and applying integration by parts to the last term in equation (*) we see that Wi, 
i E {1,2} are generalized solutions of the Cauchy problem for the generalized porous medium 
equation with non-linearity S and convection —6 in the sense of Definition 4 in IfTOl . From 
Theorem 4 in IfTOl we deduce W\ = w 2 by noting that \x and a can be easily extended to the whole 
ofR + without violating Hypothesis 1 of IfTOll . Hence, V\(t) = z/ 2 (t) foralH E [0, T] as desired. □ 
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